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Abstract 

Approximate group analysis technique, that is, the technique combining the methodology 
of group analysis and theory of small perturbations, is applied to finite-difference equations 
approximating ordinary differential equations. 

Finite-difference equations arc viewed as a system of algebraic equations with a small pa- 
rameter, introduced through the definitions of finite-difference derivatives. It is shown that 
application of the approximate invariance criterion to this algebraic system results in relations 
that can be viewed as prolongation formulae and the invariance criterion for the differential 
approximation of these finite-difference equations. 

This allows us to study the group properties of the finite-difference equations by analyzing 
the group properties of their differential approximations, which are the differential equations 
with a small parameter. 

In particular, the question of whether the group, admitted by the original differential equa- 
tion, can be corrected by adding the first-order perturbation to it, so that the resulting group 
with a small parameter is approximately admitted by the finite-difference approximation, is 
studied. It is shown by examples that, for a given differential equation, its finite-difference 
approximation and the group, such a correction may not always be possible. It is also demon- 
strated that the finite-difference approximation can be modified in such a way that the correction 
becomes possible. 



1 Introduction 

An invariance of the differential equations with respect to one-parameter groups of transforma- 
tions in the space of independent and dependent variables carries important information about the 
fundamental properties, such as conservation laws, of a physical system that these equations may 
describe. 

Knowledge of the groups admitted by the system of differential equations may also allow one to 
reduce the order of the system, to find important particular solutions, or to produce new solutions 
from a single solution which is known Q . Preserving the group properties of differential equations 
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in their discretizations seems to be a natural requirement which would ensure that the fundamental 
properties of the approximated differential equation are preserved in the discrete model. 

A step in the direction of analyzing the properties of purely difference equations was made 
by Maeda ||, where the questions of linearization and reduction of order in ordinary difference 
equations were treated. Levi and Winternitz H applied a Lie groups technique to differential- 
difference equations arising in study of lattices. Knowledge of these group properties allowed the 
authors to obtain some non-trivial invariant solutions. 

Dorodnitsyn Q developed a general theory of transformation groups in the spaces of discrete 
variables. By properly constructing the space of discrete variables and the finite-difference ana- 
logues of the differential operators, he demonstrated that the methodology results in the prolon- 
gation formulae for the finite-difference derivatives, which are similar to those formulae in the 
continuous case. This similarity allowed the author to extend some of the group-analytical meth- 
ods to the finite-difference equations. In particular, the difference analogy of the Noether theorem 
is proven. 

Unfortunately, in general, replacing the system of differential equations with their finite-difference, 
finite volume or finite element discretizations introduces noninvariance into the discretization. This 
noninvariance stems from the fact that performing the discretization requires utilizing a discrete 
set of nodes in the space of independent variables - a computational grid. Formally this implies 
that the group of transformations acting on the discrete variables needs to leave invariant not only 
the approximating finite-difference system, but also a set of algebraic equations defining the grid 
Although some of the important groups of transformations arising in applications do satisfy 
this property, generally speaking, this requirement results in the discrete system losing many of the 
groups of the original differential equation. 

The natural way to avoid this restrictiveness of the above approach is to weaken the require- 
ment of invariancy of the finite-difference system. One of the possibilities, suggested and studied 
by Shokin (7j, is to consider the differential approximation of the finite-difference system. By for- 
mally applying a Taylor series expansion to the finite-difference approximation, one can obtain a 
differential equation which can be treated as a perturbation of the original differential equation due 
to discretization. The perturbation terms, added to the original equation, depend on the step sizes 
of the grid. 

In the general theory of differential approximations of the discrete models is developed and 
group properties of the differential approximations are studied. Following this methodology one 
can find the groups admitted by the first differential approximations and deduce some important 
properties of the discrete model from analysis of these groups. 

Despite the simplification introduced by considering the first differential approximations instead 
of the system of algebraic equations, the approach 0j still requires answering some important 
questions. 

The step sizes of the grid participating in the differential approximation are the parameters 
which need to be transformed when a group of transformations acts on the independent variable. 
Since the step size is defined as the difference between the values of the independent coordinate 
at the end points of the mesh interval, it can easily be seen that the defining equations || for the 
general symmetry group of the given system will no longer be local in the sense that they will 
contain not only the derivatives of the infinitesimal operator's components, but also the values of 
these components at the neighboring nodes. Although for many important symmetry groups this 
non-locality can be avoided 0], this fact makes finding the general symmetry group as hard as 
solving the finite-difference equation in closed form. 

Another important question related to the approach 0] follows from noticing that, in general, 
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the symmetry group admitted by the differential approximation of the difference equation does 
not leave the original difference equation invariant. Therefore, further formalization is needed in 
order to study the properties of the discrete model through analyzing the groups of transformations 
admitted by the differential approximation. 

In this work the recently developed technique of approximate group analysis due to Baikov, 
Gazizov and Ibragimov Q is used as a basis for such a formalization. The discrete model is viewed 
as the system of algebraic equations with a small parameter introduced through the definitions 
of the grid step sizes and finite-difference derivatives. The groups of transformations leaving this 
system of algebraic equations together with the definitions of step sizes and the difference derivatives 
approximately invariant |2| are then sought. It is demonstrated that the groups of transformations 
satisfying these properties can also be found by considering the differential approximation of the 
discrete model, provided that the action of the group on the step sizes is extended in a special way. 
This "prolongation" of the group to the step sizes is local and therefore it is easier to solve the 
determining system analytically. 

2 Approximate symmetries 

In order to make the presentation complete and to introduce some important concepts to be used 
in the following, the basic notions of approximate symmetries are outlined in this Section following 

1- 

2.1 Approximate one— parameter groups of transformations 

Let z = (z\, zn) E \R n denote an independent variable, e be a small parameter, and p be a fixed 
natural number. 

Consider the following one-parameter family of approximate transformations acting in IR N : 

z / = f(z,e,a)+o(eP), (1) 
where f = (/i, /jy) and a is a scalar parameter of this family of transformations. 

Definition 1 We say that transformations (Q) form an approximate one-parameter group (a.o.p.g.) 
with respect to the parameter a if 

f(z,e,0) =z + ( £ P), (2) 
f (f (z, e, a), e, b) = f (z, e, a + b) + o{e p ) (3) 

and also a = if f (z, e, a) = z + o(e p ) is true for all z. 

Definition 2 An approximate infinitesimal operator corresponding to an a.o.p.g. (Q) is 

X = fd z + o{eV), (4) 

where d z = (d Zl , d ZN ) T , £ = (6, 6v) T and&(z,e) = ^/i| a =o> i = 1,...,N. 

The following analogue of the Lie theorem || for a.o.p.g. establishes a relation between an 
a.o.p.g. (||) and its approximate infinitesimal operator (||). The proof can be found in Q. 
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Theorem 1 (Lie theorem for a.o.p.g.) If the transformation (^) forms an a.o.p.g. with an 
approximate infinitesimal operator fy), then the function f(z,e, a) satisfies the equation 

^f(z J e,a)=€(f(z,e,a) J e)+o(eP). (5) 
Conversely, for arbitrary smooth function, the solution of the Cauchy problem 

z'| o=0 = z + o{e p ) 
yields an a.o.p.g. with group parameter a. 



Example ]l]], Consider the case p = 1 . Assuming that the function f (z, e, a) is analytic, write 
the a.o.p.g. (jl|) as: 



z = f (z, a) + e f (z, a) + o(e) 



(6) 



An approximate infinitesimal operator can then be written as 



X = X+eX +o(e) = £ T (z)<9 z + e£ T (z)d z + o(e) 
1 o 1 



(7) 



where 



£(z) = |-f(z,a) 
o aa o 



a=0 



€(z) = |-f(z,a) 
i oai 



a=0 



Using (^) and (^) in Theorem |] we find that f and f satisfy the equations 

l-f = *(f). 
oao oo 

^-f = £(f) + £'(f)f, 

aa l i o o o l 



3 

where = — Conversely, the functions f and f can be found from the approximate infinitesimal 

o oz o 1 

operator (|7|) by solving the following Cauchy problem: 



d 

da o 

— f 

da l 



£(f), 

o 

£(f)+£'(f)f, 

1 1 



a=0 



0. 



(8) 
0) 
(10) 



a=0 
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2.2 Approximate invariance 

Let Fi(z, e), i = 1, M be given scalar-valued functions and 

F(z,e) = (F 1 (z,e),...,F A/ (z, £ )) T . 



(11) 



Definition 3 The equation 

F(z,e)=o(e p ) (12) 
is said to be invariant with respect to an a.o.p.g. (Qj, if, for all z satisfying we have 

F(f(z,e,a) lE )=o(eP). (13) 
The following criterion of invariance has been proven in 0. 



Theorem 2 (Invariance criterion) Let the function (11) be analytic with respect to the variables 
z and e and also 



'OF 

rank [ — — (z, 0) 
1 oz 



M. 



F(z,0)=0 



(14) 



Then the approximate equation nQ ) is invariant with respect to an a.o.p.g. ([Ip with an approximate 
infinitesimal operator (^jj if and only if 



XF(z,e)\$ = o(en. 



Example Q, ||. Let us consider the case p = 1 and obtain the necessary and sufficient con- 
dition under which the equation 

F(z,e)=F(z) + eF(z)=o(e), (15) 

where F and F are both analytic functions, is invariant with respect to the a.o.p.g. (||), (|^). First 

notice that substituting (||) into (p^), expanding the result in powers of e and taking into account 
analyticity of the involved functions results in the following set of "exact" equations, which is 
equivalent to the invariance condition (|ll|) with p = 1: 

F(f(z,a)) = 0, 
F(f(z,o)) + F / (f(z,a))f(z,o) = 0, 

where F = tt-F. 
o oz o 

Suppose that F satisfies the nondegeneracy condition following from ( |14D . Then the infinitesimal 
criterion of invariance can be written, according to Theorem B as follows: 

= o(e), (16) 

(15) 



X F(z) +eF(z; 

,0 1 
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where X is given by (^). Following Jp, [g], we substitute z = z+ez + o(e) into ( |l6| ) in order to 
obtain the equivalent set of exact equalities: 

XF(z) = 0, (17) 
o o o' 

XF(z) +XF(z) +z9 z (XF(z)) = 0. (18) 
l o o o i o i o o V 

The above two equalities must be satisfied when 

F(z) = 0, (19) 
F(z) + zd z F(z) = 0, (20) 

as follows from (|i~5l). 

In the following we further restrict our consideration to the case M = 1, that is F = iq(z,e). 
In practice, as suggested in Q, it is convenient to get rid of z in the above relations by 

employing the following fact. It can easily be seen from 

(0) and © that X is an infinitesimal 

operator corresponding to the exact transformation group which leaves equation ( |l7|) invariant. 
Furthermore, as follows from Proposition 2.10 from p. 84], the left-hand side of (|l7]) can be 
written as follows: 

X F(z) = A(z)F(z), (21) 
o o o ; V o (r 



where A(z ) is some smooth function. Then, using pQ ) and substitutin g (pi|) into fll8|), we rewrite 
the determining system (|l7|), ( |i~8|) together with the relations ( |l~9|) and (|20|) as follows: 

XF(z) = A(z)F(z), (22) 



XF(z)+XF(z) = A(z)F(z). (23) 

10 1 1 

If the invariance group for the ("non-perturbed") equation F = is known, the relation ([22]) can 

be used to determine the smooth function A(z). This function is then substituted into (^), which 

must be satisfied by virtue of (fjl^). Solving (p3|) together with (|i~9|), one finds the "perturbation" 

X of the infinitesimal operator, 
l 



Example H. Letting p = 1, N = 2, z = (z\, z-i) , equation (15) reads 

F(z,e) = F(z) +eF(z) = z\- l + e(z%lnz 2 - z\) = o(e), 
o 1 

The above equation can be considered as a small perturbation of the equation F(z) = 0. The latter 

o 

equation obviously admits the group X = d Zl . One may consider a small perturbation of this group 

given by a.o.p.g. X = X+£X+o(e) and establish the condition on X ensuring that the a.o.p.g. 

leaves the equation invariant. 

First notice that (^) results in A(z) = 0. Using it in (^), one arrives to the condition 



i 

which needs to be satisfied when (z 2 ) : 



(6(z)) 2 = (zi) 2 , 
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One a.o.p.g. satisfying this condition is X = (1 + ezf)d Zl + eziz 2 d Z2 +o(e). Solving the Cauchy 
problem (§)-(|lOl), one can find the transformations corresponding to this a.o.p.g.: 

a 3 

z[ = Z\ + a + e{az\ + a 2 z\ + — ) + o(e), 

o 

2 

z 2 = z 2 + e(az 1 z 2 + — — ) + o(e). 

3 Approximate symmetries of finite— difference equations 

3.1 Finite-difference approximation as a system of equations with a small pa- 
rameter 

Consider an ordinary differential equation of order n: 

F(x,yj 1 \...,y^) = 0, (24) 
y = y(x), x e (a,b). 

Introduce the grid as the set of points x a € (a,b), a = 1, N, and the grid function u a , a = 1, N 
defined at these points and to be used to approximate the values of y = y{x) at the nodes of the 
grid. 

Consider also the set of algebraic equations with respect to the unknowns x a and u a : 

L a (xi, ...,x N ,ui, ...,u N ) = 0, a = l,...,N. (25) 



In the cases of interest to our work, this set of equations approximates the differential equation ( p4[ ) 
in the sense defined below. 

In order to introduce a small parameter into (p5|), extend the system of equations ( p5| ) by 
introducing the new unknowns and adding the following equations to it: 



X a +1 


x a 


= eh a 


(J) 




p 




= E 

k=l 




n (0) 


= Ua', 


N-l, 


1 

a 


= I,- 



e 

ft!' 



(26) 
(27) 

(28) 



a = 1, N — 1, a = 1,...,N, p > n, J = 0,...,p, 

where e is a small parameter and the added unknowns, h a and Ua \ (J = 0, 2p) will be referred 
to as, respectively, step sizes and finite-difference derivatives of J-th order. 

Although this system is comprised of algebraic equations, equations ( |27| ) can be viewed as a 
result of a formal expansion in the Taylor series at the points x = x a , which is used in order to 
obtain the values of where (3 = 1. The following two Lemmas demonstrate that, by virtue of 



the equations (26), ( [27D and (^), this analogy with the Taylor expansion can be extended to the 
cases when x a and x a+ p are two arbitrary nodes of the grid, and that it can be extended to the 
smooth functions of the variables x a and u a . 



Lemma 1 If 1 < a < N , (3^0, l<a + /3<NandO<J<p-l, then, by virtue of (\2Q), (g) 
and 

u% - <4 J) = E i h k a+ ^ a u^ + o(sV- J ), (29) 
k=i K - 
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where 



( a+P-l 

E h i iff 5 ^ 1 



h, 



a+f3,a 



i=a 
a-l 



E h i ifP<-i 

i=a+f3 



Proof. Consider first the case /3 > 1 and use induction by (3. If (3 = 1, 

(ID- 



is true because of 



Assume (|29|) is true for some (3 = I > 1. Then, using (|26|), (|27|) and (|g), we get by direct 
calculation (see Appendix [A| for details): 



(J) 

a+l+l 



U 



(J) + faO _ J J ) ) i f n ( J ) _ u 0/h 



(J) + E^' a+w /i J+,) +«r J ) 



(30) 



which proves the validity of (|29| ) for /3 = / + 1. Therefore, by induction, the Lemma is true for 
/3> 1. 

Consider now the case /3 < — 1. We again perform induction by (3 starting from (3 = —1 and 
going in the direction of decreasing (3. To prove the statements on each of the steps of this induction, 
induction by J will be used. 

To prove (09) for f3 = — 1, consider induction by J from J = p to J = 0. First, for J = p 
(|27|) implies = + o(l), which means that ([29|) is valid when /3 = — 1 and J = p. Assume 
that (||) is true for (3 = -I and J = J + 1, Prove it for J = J . Using (|(|, (§7|), (||) and 
induction hypothesis one obtains by direct calculation (see Appendix |A] for details): 



u 



(Jo) 
a-l 



,(Jo) 



,(Jo) 



P-Jo e k h k 



E 

k=l 



a-l (Jo+A:) 



A;! 



P-Jo -kuk 

e n a _ 1 



E 

k=l 
p-Jo 



hi 



<-{' + o(e p - Jo ) 



m=0 



+ o e 



(31) 



4 Jo) + E ^^f^^ Jo+i) 



+ o(e p - Jo ), 



and therefore (|29| ) is true for = — 1 and all possible values of J. 

Now, assume ( p9| ) is true for all f3 > 0q, where (3q < —1. Prove it for j3 = (3^ — 1. Again 
use induction by J. For J = p (|29|) immediately follows from (|27|). Suppose ( p9|) is valid for 
J = Jo + 1, ■~,p. Using (|26|), (p^), (p8|) and induction hypothesis prove that it is true for J = Jq 
(see Appendix [A|) : 



u 



Jo 

a+0 o -l 



Jo e fc (-/i Q+/3o _i) fc (Jo+fc) 



E 

fc=0 



fc! 



P- J o fe/ h \k 
^_y~ n a+f3 ~l) 



k=0 



k\ 



p—Jo—k c-mu 

E 

m=0 



a +Po, a „.(Jo+k+m) 

i "a 
m! 



+ o(e p - Jo ) 



(32) 



,(Jo 



p-Jo 



1=1 



' 0) + E 7F^ Jo+4) ^o-i, a + °(^" Jo )' 
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which proves fl29|) for all possible /3 < — 1 and J. 



□ 



Lemma 2 If 1 < a < N, and £(x, u) is a p times continuously differentiable function, then, by 

. -D m Z(x a ,u a ) + o(eP), 



mrfue o/ (|F|) emd (gj 

£(x Q +i,u Q +i) - £,(x a ,u a ) = E 



m=l 



ml 



(33) 



where D = d x + u^c^ + u^d u (i) + 



Proof. Using smoothness of the function £(x, u) and rearranging terms in a Taylor expansion (see 
Appendix [B|) write 



m=0 



where 



???.! 



m /to— l / -\j \ 



i=0 \n=0 



+ o(e p ), (34) 



(0 

(n) 



h+...+l„=l 



U = u®/l\ and L?i = <5 Z)0 . To prove the Lemma, show that the expression in square brackets in 



( 34|) is equal to 



D m i(x a ,u a ) 



d x + Yu {J+1) d u 



(J) 



i=0 



E- (J+1) ^ 



(j) 



<9*£(x a ,u«)(35) 



Comparing equations (35) and (|34|), one can see that it is sufficient to prove that the expressions 
in round brackets in these equations are equal for all values of I = m — i varying between and m, 
ic 



n=0 V J / 



(36) 



This equality is proven using induction by I. When 1 = 0, the equality becomes an identity. Assume 
that the equality is true for I = k < m. By applying the differential operator J2 j u J ^ U ( J ~> to 
both sides of ( |36| ) written for I = k and manipulating the terms in the left-hand side of the resulting 
equation (see Appendix |b|), one obtains equation (]36| ) with I = k + 1. This completes the proof of 
the Lemma. 
□ 
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Consider one of the equations ( |25|) corresponding to a = ao- Using equations (|26|), ( p7| ) and 
Lemma |] equation (^5|) may be rewritten as follows: 

L ao (xi, ...,xn,ui, ... ,un) 

(a -l jV-1 
X a Q S ^ h{, . .., OCcto j ~t~ ^ ^ ^ ^i; 

i=l i=ao 

P k "0-1 V k N-l \ 

^o + E at- E M fe * •••> «oo. - ^o + E^[(E +°( £P ) ( 37 ) 

k=l ' i=l k=l ' i=ao J 

= 4° )(x ao ,<),u«,...,«g,/ il ,... ! / li v-i)+E^ 

where 1 < r < p is some integer number. 

Definition 4 T/ie system of algebraic equations (pffi), (U^ and (2<5j is said to approximate 



equation (HJj a£ £/ze node x ao , if, for all values of x aQ , u^J, Ua[~^..., u$}, hi, h^-i, we get 
L {0) (x u (0) u {p) hi hATi)-F(x u (0) u W u {n) ) 

^ao K^acn a a ' "ao ' •"' "ao ' lj '"' "'N-l) — V-^ao' "ao ' "ao ' ""' "ao /' 

where F is the function from (f^J. 
Definition 5 A differential equation 

F(,,y,^,... J W) + ^ ! 41(^l/ I 2/ W 1 .-- 1 !/ (p Ui,.,^-i) = 0, (38) 

i=l 

where 1 < r < p, e is a small parameter and hi, /tj\r-i are parameters (not small) is called a 
differential approximation ofr-th order at the node x ao for the system of equations $2$), $2(\), W\ ) 
and $2t\). If r = \, the differential equation $3$ ) is called the first differential approximation at the 
node x ao . 



Assuming that the set of equations (f26[), (27) and (28) approximates the differential 



equation (p4|), one can rewrite the former set of algebraic equations as follows: 

f{Xa>,U a , ,U a , ,...,U a , ) + 2_^ £ L a> \ x a' ,U a , , U q/ ,...,U a , J 



fc! 



k=l 



(0) 



(O, 


(39) 




(40) 


o(e r ), 


(41) 


0, 


(42) 



where a = 1, N — 1, a 7 = 1, ...,N, 1 < r < p, p>n, J = 0,...,p. 

As a result of the above, we obtain the set of equations (|39|), (|40|), ([4l|), (p^ ) with a small 
parameter e. In the following, approximate symmetries of this system and their relation to the 
symmetries of the differential equation ( |24| ) will be studied. 
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Example. Consider a general first order differential equation resolved with respect to the first 
derivative: 

F(x,y,y W ) = y W - f(x,y) = 0, x £ (0, 1) (43) 
and its finite-difference approximation based upon forward differencing: 

L a (x a , x a+ i, u a , u a +i) = Ua+1 _ U ° - f(x a ,u a ) = 0, a = 1, ...,N - 1. (44) 

Xa+1 X a 



Set p = 2 and obtain the system (|39|), ([10|), (41), (|42j ) corresponding to the finite-difference equation 



fl4§. Using (gg) and @ in @, we get: 



^ 1) -/(x a ,n(°)) + ^n( 2 ) = o(e), 
x a +i-x a -eh a = o(e), 

Q + l 



,(°) 



0, 



where a = 1, ...,JV — 1, a' = 1, N, J = 0,1, 2. 

According to the definitions given, the above system of finite-difference equations approximates 
a differential equation (^) at all nodes x a , a = 1, N — 1. At any of these nodes the first 
differential approximation corresponding to (E3) is 



y (1) -f(x,y) + £ -^y {2) = o. 



As can be seen from this Example, the term with a small parameter, appearing in the first differ- 
ential approximation, is linear with respect to h a . The following Theorem gives the necessary and 
sufficient condition that equation ( |38|) with r = 1 needs to satisfy in order to be a first differential 
approximation of some finite-difference approximation of equation (|24|), 

Theorem 3 Let N > p > n, then the equation 

F(x,y,y^\...,y^) + eL^(x,y,y^\...,y^,h 1 ,...,h N . 1 ) = 0, (45) 

is a first differential approximation of some finite-difference approximation of equation at the 
node x ao if and only if 

X(eLW) = 0, (46) 

where 

n 

X = -ed E + ^h CH d CH . (47) 

i=l 

Proof. 1) Let the condition ([!(]) be satisfied. The following describes the procedure resulting in 
a finite-difference approximation which has ( |45| ) as its first differential approximation. 

The differential operator X in ( f47|) can be treated as an infinitesimal operator of the o.p.g. 
acting in the space of p + N + 2 variables e, x, y, y^\ y^ p \ hi, /tjv-i- As one can easily see, 
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x, y, yM, ...,y( p \ehi, e/tjv-i are p + N + 1 functionally independent invariants of this o.p.g., and 
eLao is also an invariant of the same group because of fl46|) . Using Theorem 2.17 (||, p. 88) and 
formally substituting x a+ \ — x a for each expression eh a and x ao for x, we establish that 

eL a] = h{x,y,y^ ] ,-,y (j '\ehi,...,eh N ^i) = f 1 (x ao ,y,y {1 \...,y^ p \x 2 -xt,...,x N - x^-i) 

= /i(y,y (1) ,---,y (p) ,^i,---,£A0 

for some smooth functions / and /. 

Utilizing ( |29] ) and the fact that N > p > n substitute y, y^, y( p ) with their respective finite- 
difference approximations of the first order of accuracy on the grid composed of the nodes x%, xn 
with u a denoting the value of the grid function at the node x a : 

-,y ip) -,xn) = fi(u 1 ,...,u N ,x 1 ,...,x N ) + 0(1) 

for some function f\. Acting in a similar way, substitute the arguments y, y^\ y^ n ' with 
their finite-difference approximations of the order higher than the first in the expression for 
F(x,y,y^ , ...,y^), and also substitute x with x ao in the same expression: 

F(x,y,y W ,-,y {n) ) = fo(u 1 ,...,u N ,x 1 ,...,x N ) + o(e) 

for some function /o. 

Now chose the finite-difference approximation at the node x aQ to be 

L(m, ...,u N ,xt, ...,x N ) = f Q (ui, ...,u N ,xi, ...,x N ) + A(ui, ...,u N ,x 1 , ...,x N ) = 0, (48) 

where fo and f\ are the previously constructed functions. 

Reversing the above speculation, we can see that, by virtue of (29) and (26): 

L(ui,...,u N ,x 1 ,...,x N ) = F(x ao ,u ao ,u^, -,u^J) + eL^(x ao ,u ao ,u^,...,u^,hi,...,h N ^ 1 ) +o(e), 



that is, ( f45|) is a first differential approximation of 

2) Let (|25D with a = ao approximate (p4j) . Introduce the following notation: 



L ao ( x *o > «£! i «§! ,..;uj$,eh 1 ,...,eh N - 1 ) 

(ao-l N-l 

X ao £ ^ ] Xqq, 3?q.q -|- £ ^ ' hi, 

i=l i=ao 

V £ k "0-1 P £ fc W-l 

u *o + E tj(- E k) fc «2?> ■••> u °o> + E jfef (E h *) fc "2? 

fc=l ' i=l fc=l ' i=ao 

Notice that it follows from this definition and (|37l), that XL ao = and also that La = (0eA*o)le=O' 
Using these relations, write 

n n 

E h a>id ai (d E L ao ) = 9 e (E h ai d ai L ao ) = d £ (ed £ L ao ) = d £ L ao + ed £ L ao , 
i=l i=i 



Taking the limit e — ► in the latter equation results in (47). 
□ 
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3.2 Approximate groups leaving finite— difference approximation invariant 

Consider a one-parameter group of transformations acting on the variables x and y and represented 
by its infinitesimal operator X = t;(x,y)d x + ip(x,y)d y . Action of this group is extended to the 
derivatives yW, ...,y^ using prolongation formulae Pj: 



X = £0, y)d x + <p(x, y)d y + V (i) d y 



(0> 



i=l 



where 



•,(0 



^" = Dtp*-'- 1 ') -y®D£, 
D = d x + y^d y + y^d yW +..., 
(o) 



(49) 



(50) 



ip. 



Action of this group of transformations can also be naturally extended to the "finite-difference" 
variables x a and u a , (a = 1, N) as 



N 



(51) 



a=l 



In order to extend action of this group to the variables , ...,u^\ (J = 1, 2p) and hi,..., h^-i 
while preserving "locality" of the prolongation formulae, consider an a.o.p.g. represented by its 
infinitesimal operator 



X (A) = ^ £ fc x (A) +o(£ r )! 
k 



(52) 



X 

k 



(A) 



fc=0 
N 

E 

a=l 



2r 



£(x a ,u a )d Xa +tp(x a , u a )d Ua + V tptyd u) 

7. 7_ ' " 7 "TV 



8=1 



AT-1 



+ E ?A- 



a=l 



In the following, we will be interested in those of the groups (f)3), which leave invariant the set 
of approximating finite-difference equations (|39|) , ( |40| ) and ( |4~l| ) . The following Theorem establishes 
the analogues of the prolongation formulae ([3(]) for the finite-difference equations: 



Theorem 4 Equations j\4u[ ) and (41) are invariant with respect to the a.o.p.g. (5i) if and only if 

(53) 
(54) 



a a 
k 



k+1 U 

1=1 l] *+!-! 



i J+1) = o^-t4 J+1) ^(*«,««), 

fc fc 



where < k < r — 1. 

Proof. Using the infinitesimal criterion of approximate invariance (Theorem 0), apply the operator 



p2|) to equation (40) to obtain: 

r 

£ a+l - £ a + E 6 (£ "+ 1 ~~ ^ 1 



fc=l 



fc-1 



o e 
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Using Lemma 0, rewrite the left-hand side of the above equation as follows: 



E z k % «+i - i «) - E c k k \ a = E E f«-E e Vi a+0 

fc=o fc fe fc=i fc=o/=i l - k k=i k 1 



>(0 



4 »J 



1=1 u=l 



- * a - O a } + o(£ 
11 i_ Z i-l 



which immediately results in (p3[). 

Similarly, having applied Theorem ^ to equation (f4l|), one obtains 



fe=0 



1=1 



i-l 



(i - 1)! fc 



(J+i) 



which must hold by virtue of (^) and . Applying Lemma Q and ( |53| ) , transform the left-hand 
side of the above equation as follows: 

M! t a i\ t a {i-l)\k aUa \ 



i=l k=0 
r I 



h 



i-l 



E^|^-f tf« (i-l)! A 



+ o e T 



i=l i=l 
r I 



h 



;=i i=i 



»i i 



l-i 



n=l 



Equating the obtained expression to o(e r ) results in the following set of equalities: 



„-_1 | l—i l — i 11 = 1 ^ — * 



0, 



i=l 



(55) 



where I = 1, 2, r. 

Use induction by 2 to prov e (|54] ) for all < k < r — 1. Setting / = 1 in fl5"5|) immediately results 
in (|54| ) with fc = 0. Assume (|54|) holds for all k = 0, 1, fci and prove it for k = k\-\-l<r — 1. 
Consider (pq ) with I = k\ + 2 . It is easy to see that the terms in the brackets, corresponding 
to i = 2,...,ki + 2 nullify because of the induction hypothesis, while the remaining term, being 
equated to 0, gives (54) with k = k± + 1. 
□ 



One can immediately notice the similarity between the prolongation formulae for the finite- 
difference derivatives (|54]) and their continuous analogues (|50|). This similarity allows one to an- 
alyze and to establish some important group properties of the finite-difference approximations in 
terms of the group properties of their differential approximations. Really, consider also the sys- 
tem of algebraic equations with a small parameter (|39l), (]40|),(|4l|) and ( |42"1) resulting from some 
discretization of the equation ([24]) . According to Theorem the necessary and sufficient condition 
for the subsystem (^),(^Tj) to remain invariant under the action of the group ([52]) is that ip^ and 

k 

a a are given by ([53]) and (|5^) . Using these expressions and applying the operator (|^) to equation 
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) , one can see that the obtained system of determining equations is the same as the one resulting 
from the application of the infinitesimal criterion of invariance to the differential approximation 
provided that action of the group is extended on h a according to fl53|). 



Example. Consider again an ordinary differential equation (|2J) and the approximating set of 
algebraic equations (j39|)-(|42|) with r = 1. The first differential approximations corresponding to 
the nodes x a , a = 1, N read: 



F(x,yJ 1 \...J^) + sL^(x,yJ 1 \...,y^,h 1 ,...,h N ^) 
Consider the following a.o.p.g. 



0. 



X 



where 



X 

k 

(0 



1 / N-l \ 
k=0 \ k *=i k J 



2p 

£(x, y)d x + ip(x, y)d y + J2 <P W ^w > 



(56) 



(57) 



k k 

fc+i h i 



(58) 



k 



1=1 



k+l-l 



that is, in addition to conventional prolongation formulae (p8|), the relation (53) is used to extend 
the action of the group on the parameters h a . 

Using Theorem |2| and the methodology outlined in the Example following this theorem, write 
the necessary and sufficient conditions, under which equation (|56|) is invariant with respect to the 
a.o.p.g. ([57]), in the following way: 



XF-XF 
o 



N-l 



XF+(X+Y,Z°dhjLW-\LV 



a=l 



o. 

0. 



(59) 
(60) 



F=0 



where A = A(x, y, y^ 1 ' , y^) is some smooth function and the fact that F does not depend on h a 
has been taken into account. It follows from Theorem ||] that L^ 1 must be a linear function of the 
parameters h a . Using this fact and also noticing that the first term in ( |60| ) does not depend on h a 
and therefore should be equated to zero separately, rewrite (59) and ( p0| ) as follows: 

0, 
0. 



XF-XF 
o 



XLU + {D£-\)lW 
o o 



F=0 

(XF)|f=o 



0, 



(61) 

(62) 
(63) 



The condition ( pT| ) implies that X must be admitted exactly by the original differential equation 
(24). This condition also allows one to determine the function A. If the additional condition (p2|) 
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is satisfied, the group X can be corrected by adding a small perturbation eX to it, so that the 

o 1 

resulting corrected group leaves the first differential equation (^6|) invariant. Finally, ( |63|) shows 
that the correction X may be any group admitted by the original differential equation (pif), or it 

can also be an operator with all zero components. 



Definition 6 A group fclQ), leaving the differential equation invariant, is called stable with 
respect to the approximation resulting in ftjfyj ( EjP an d §4^) if and only if there exists an in- 

finitesimal operator X, such that the corresponding a.o.p.g. X^ + eX' A ' + o(e) leaves the system 
l ' 1 

of equations ftsQ), $4~(\), and ffify invariant up to the o(e) terms, ie the a.o.p.g. x' A ' + eX ^ + o{e) 



leaves invariant the system (pty, (4L ), (j4A) an d (4^)> * n which r is set to be 1. 

Using this notion, the results of this Section can be summarized in the following Theorem. 



Theorem 5 A group leaving the differential equation (24) invariant is stable with respect to 
the approximation which has the first differential approximation (ffli) if and only if 



xl$ + (d(Z)-\)lW) 

o o / 



F(x,y,yM ,...,y( n ))=Q 



(64) 



for all 1 < a < N . 



3.3 Examples 

Example. Consider the following ordinary differential equation: 

F(x, y, y') = {y- x)y' + x + y = 0, 
y = y{x), xe(a,b) 



and the group 



X = yd x - xd y - (1 - (yT)dy, - 2,y'y"d y 



Since XF = —y'F, the above group leaves the equation invariant and A = —y'. 

Consider also the following finite-difference approximation of this differential equation: 

{y a ~ x a ) Va+1 - Va + x a + y a = 0, 

Xa+l X a 

where x a € (a, b) and a = 1, N. 

The first differential approximation corresponding to this finite-difference approximation reads: 

(y - x)y' + x + y + e%y"(y - x) = o(e), 



and therefore, L« = ^fy"{y — x) 
From the above we get: 



hr, 



XLW + {D (0 - A)L« = -'-fy"((y - x)y' + x + y) = 0, 
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if (y — x)y' + x + y = 0. 

Therefore the considered group is stable. 



Example. The ordinary differential equation 

F(x, y, y') = y' - (x + l)y 2 - y = 0, 
y = y(x), xe(a,b) 

is invariant with respect to the group 

X = (y + xy 2 )d y + (y' + y 2 + 2xyy')d y > + (y" + 4yy' + 2x(y') 2 + 2xyy")d y », 

because XF = (1 + 2xy)F, and therefore, A = 1 + 2xy. 

Consider also the following finite-difference approximation of this differential equation: 

{x a + l)y 2 a -y a = 0, 



where x a € (a, b) and a = 1, N. 

Applying the procedure described above, one obtains the following first differential approxima- 
tion corresponding to this finite-difference approximation: 

y' -{x + l)y 2 -y + e-^-y" = o(e), 

and L { a } = \y"{y-x). 
From the above we get: 

XL« + (D(0 - A)L« =h 2 a (l + (x+ l)y)(2 + x + xy{x + 1)), 

and the expression on the right-hand side is not identically zero when y' — (x + l)y 2 — y = 0. 
Therefore the considered group is not stable. 



Example. Build a finite-difference approximation for the differential equation in the previous 
example, such that the group considered in the same example becomes stable. 

First, find a first differential approximation of such a finite-difference approximation: 

y'-{x+ l)y 2 -y + eL« (x, y, y' , y", h a ) = o(e) 
La can be sought as 

L { aHx,y,y',y",h a ) = -yy" + h a L*(x,y). 
The stability condition ( |64| ) then results in the following equation for L*(x,y): 

(1 + xyf ( y(1 + a . y) ) +(! + (* + l)y)(2 + x + xy(x + 1)) = 0, 

which gives 

L* = -^±^V(1 + xy) - ^ + C(x)y(l + xy), 
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where C(x) is an arbitrary function. 

Therefore, the first differential approximation of the stable finite-difference approximation must 
read: 

y'-(x + l)y 2 -y + ehd^- (x + 1} V (l + xy)-^ + C(x)y(l + xy^j = o(e). 

In the previous example we observed that the term with the second derivative in brackets ap- 
pears as a result of approximation of the first derivative using forward differencing. This observation 
leads us to the following finite-difference approximation of the given equation, which has the above 
first differential approximation: 

Ua+l ~ ya / , i\ 2 | / \ ( ( x a + I) 2 2 /-i , \ Va , m \ n , \ 

(Xa + l)y a ~y a + [xa+i - x a ) y a (l + x a y a ) j + C(x a )y a (l + x a y a ) 

x a+l x a V x a x a 

As the above consideration demonstrates, adding the last term to the finite-difference approx- 
imation results in the considered group becoming stable. 

A Calculations from the proof of Lemma 1 



Derivation of equation (p0|): 

,(J) -„xj)m(J j I _„( j )u^ j )._„w 



a+l+l — U a + {U a+l+l U a+ i) "I" "a ) 



P e k h k 



i J) + E 
fc=i 



Q + Z (J+fc) 



P-J £ k h k 



k\ ~ Ua+l ' H. fc! 



+ E 



fc=l K - 

P-J Jc r (p-J-k mum 

h k I (J+fc+m) I , > fc (J+fc) 



(J) + E 
fc=i 



m=0 



??).! 



+ o(e 



p—Jp—J—k k+m P~J k 

Zv Zv fc!m! n a+l n a+l,a u a T Zv fc! n a+l,a u a T °l £ J 

fc=l m=0 ' ' fc=0 

p—Jp—J i p—J i 

Z^ Z^ _ ^ «! a + l > a a + ^ j 

fc=l i=fc ' i=0 

p-Jp-J a 

Z, Z, fc!(i _ fc)! '»a +I 'Wj,a«« + °^ ) 
fc=0 i=k v 



To prove ([H]) use ([26|), (p^) and the induction hypothesis: 

P-Jo c kuk 

E% i ^ fc) + ^ p - Jo ) 



(Jo) 



a-1 



U 



(Jo) 



k=l 
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P-Jo c kuk 
k=l 



P y° e m (-h a ^ 



-u 



(Jo+k+m) 



+ o(e p - Jo ) 



P-Jop-Jo u(k+m) 

= ^ Jo) -E E e^H-ir^i-u^+o^) 

k=l m=0 



k=l m=0 
p~J p-J ,j 

= ^ Jo) - E E ^-iy- k j^^^ + 

= U W _ £ e i(_i)ii^2=l u (Jo-M) 

i=k 



il 



Y(-i\ k — _ — 

iti k\(i-k)\ 



+ o(e- 



.p-Jo^ 



= ^ Jo) - P J2 ^(-^^r^a ^ [(i - i) i - 1 

i=A; 



+ o(e 



.p-^o> 



= ^ Jo) + E° ^ ^' l)t ni Jo+t) + o(e p - Jo ), 
Similarly, utilizing fl26|), (27), ( |28| ) and the induction hypothesis, we derive (|32"|): 



,M>) \ - g (-ha+00-l) (Jo + fc) -/ r p-J N 

k=o K - 

V~Jq Jzt l \k 

E g (-ft a +/3 -lj 
A;' 

k=0 

— J p-J k 



p—Jo—k mum 

y> n a+/3o,a u (J 0+ k+rn) 

m=0 



ml 



+ o(e p - Jo ) 



jfc=0 i=fc 



= E° ^ Jo+i) [E 

i=o Lfc=o ^ ' 



i=0 

P-Jo _i 



+ o(e 



P-Jo i 

= E - u a° +i) (- h *+P 
i=0 ? - 

P-^o i 

= ^ Jo) + E U J 
.•_ 1 



_i + h a+l 3 , a y + o(e p 
.p-J<n 



B Calculations from the proof of Lemma 2 

To obtain equation (j34|), use the notations introduced in the proof of Lemma |2| and write: 



E(^)%) =EW i w+°(^ 



(n) 

\Z=1 / Z=n 

Using this relation and changing the order of summation, write 



£(x a+1 ,u a+1 ) = EEwi^i 



1 J^fc^f^ajUa) 



fc=0 i=0 



kl I * dx l du k 



(x a+ i - x a y(u a+1 - u a ) k 1 \ + o(e p ) 
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k=Oi=0 { \l=l / 



= EE E h c ^% a ^ {eK)l+lL{ l^ + o{£P) 

k=0 i=0 l=k—i 

E (eh a ) m J m\ . d k £(xg, u a ) (m _ 1 . 

„ m! l^f- k\ k datdu*-* L (*-t) ( +0{£ > 
\ i=0 k=i ) 

m /m—i \ 

E [ E Tj ^(n) ^" J ^£(^0:, 
\n=0 / 



m=0 

' m! 

m=0 



+ o(fi*). 



=0 U = \n = / 

Prove now equation ( |36| ) with / = A; + 1 using the fact that the same equation is valid for I = k. 
First notice that, as follows from (|65|), I/fK can be expressed as follows: 

\ \l=l / / a =0 

Having applied the operator J2j U ^ J+1 ^ U ( J ) to both sides of equation ( |36| ) with I = k, use ( |66|) and 
the expressions for L) k ( and L^ j to transform the left-hand side as follows (the arguments of £ are 
omitted for brevity): 



J U=0 J n=0 



n=0 V J / 

= E S I 4*-d™ (1) + E « (J+1) aK + 
,,=i l / ) 

X Mn d k \ f^- 

n=l 
n=l 

= E 



n=l ^ \i=l / »=0 J a=u 

_ * *in^ iff ^v^f ^L u( a ^ 



^ +(u (i) ) fc+i a «+^ 

n=l ^ \»=1 / l= ^ ' J a=Q 

i=0 



, n! 

n=l 
A- 



(fc + l)!^ 1 \ ^ »! / J 

n=l n=0 

and therefore, equation (36) with i = k + 1 



holds. 
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